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q ■ 1 INTRODUCTION 

■ In our quest to solve the mystery of the observ ed cosmic accel- 
' eration dRiess et al.lll998l ; iPerlmutter et al.lll999D . galaxy cluster- 
ed ■ ing plays an increasingly important role as a probe of both dark 
energy and gravity, the two main classes of possible explana- 
tions for cosmic acceleration. Current data from t he Sloan Digital 
Sky S urvey (SP SS) Data Release Seven (DR7) jAbazaiian et alJ 
2009), WiggleZ telake et alJ 120091) . and BOSS feisenstein et"ai] 
201 If) are allowing us to place very useful constraints on dark en- 
ergy. The planned space mission Euclid^ will survey ~ 60 million 
emissi on-line galaxies at 0.7 < z < 2 over 15,000 square de - 
grees dCimatti et al]|2009l ; IWang et alj|2010l ; lLaureiis et ai]|201lh . 
and provide potentially revolutionary bounds on the nature of cos- 
mic acceleration. 

The SDSS data have be en analyzed using both the power 
spectrum method (s ee, e. g., iTe gmark et a l.l |2004| ; iHutsi 2005; 
Padmanabhan etal] 120071; iBlake et alJ 120071: IPercival etail 
20071 l20ld: iReid et al.l [2010l: iMontesano et al.l l201lh. and 



ABSTRACT 

We present a simple and efficient phenomenological model for the two-dimensional two- 
point galaxy correlation function that works well over a wide range of scales, from large scales 
down to scales as small as 15 h~ 1 Mpc. Our model incorporates nonlinear effects, a scale- 
dependent galaxy bias on small scales, and allows the redshift-space distortions to be scale 
and direction dependent. We validate our model using LasDamas mock catalogs, and apply it 
to the Sloan Digital Sky Survey (SDSS) DR7 Luminous Red Galaxies (LRGs). Using only the 
monopole and quadrupole of the correlation function measured from the SDSS DR7 LRGs, 
we obtain significantly improved measurements on H(z), Da(z), and (3(z): H(z)r s (zd)/c = 
0.0423 ± 0.0027, D A (z)/r s (z d ) = 6.77±0.17, and (3(z) = 0.348 ±0.058 at z = 0.35, using 
the scale range of 25 < s < 120/i _1 Mpc. We expect our model to be useful in tightening dark 
energy and gravity constraints from the full analysis of current and future galaxy clustering 
data. 

Key words: cosmology: observations, distance scale, large-scale structure of Universe 
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I Chuang. Wang. & Hemanthi] 1201 2t ISamushia et all 1201 ll; 
iPadmanabhan et alj2012f) . Although these two methods are simple 
Fourier transforms of one another, the analysis processes are 
quite different and the results cannot be converted using Fourier 
transform directly because of the finite size of the survey volume. 

The power of galaxy clustering as a dark energy probe lies 
in the fact that the Hubble parameter, H(z), the angular diam- 
eter distance, Da(z), can in principle be extracted simultane- 
ously from data through the measurement of the baryon acous- 
tic oscillation (BAO) scale in the radial and transverse direction s 
jBlake & Glazebrookll2003l ; ISeo & Eisens"teinll2003l : IWang||2006l) . 
The inclusion of information from full galaxy clustering goes be- 
yond BAO only, and enables significantly enhanced constraints on 
H(z) and Da{z). Most importantly, it allows the measurement of 
the growth rate of cosmic large scale structure, f B (z) = f3(z)b(z) 
(w here j3(z) den otes the linear redshift-space distortion (RSD) fac- 
tor l l Raised 19871) . and b(z) denot es galaxy bias), required for usin g 
galaxy clustering to test gravity dGuzzo et al.ll20 08; Wang 2008). 
While it will be possible to me asure f q (z)a g(z) dSong & Percivall 
120091) or f a (z)cr m {z)/r a {z d ) A jWandl2012l) along with H(z) and 
Da(z) in a model-independent manner from future data, current 
data require assumptions (e.g., CMB priors) being made for such 
measurements. Measuring H(z), Da(z), and /3(z) is comparably 
easier and more straightforward, and provides a stepping stone for 
future studies. 
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In IChuang & Wang] feoill) . we made significant improve- 
men ts in modeling ga l axy clustering from previous stud- 
ies ( Okumura et ai] l2008l ; ICabre & Gazianagal 120091 : | Kazin et al.l 
l2010bl) , and succeeded in making the first simultaneous measure- 
ments of H(z) and Da(z) from data, usin g the full 2D correlatio n 
function of a sample of SDSS DR7 LRGs ( Eisenstein et al.l 2001). 
and with out assuming a dar k energy model or a flat Universe. Very 
recently. iReid et ail d2012l) measured H(z), Da(z), and growth 
constraints at z = 0.57 from the monopole and quadrupole of 
the 2D 2PCF of the SDSS III Baryon Oscillation Spectroscopic 
Survey (BOSS) (SP SS DR9) sample of galaxies, assuming CMB 
priors. Most recentlv. lXu et alj d2012l) measured H(z) and Da{z) 
at z = 0.35 from the SDSS DR7 LRGs by applying density-field 
reconst ruction to an anisotropic analysis of the BAO peak. 

In lChuang & Wand d2012l) . we extended our method by ex- 
ploring the use of the multipoles of the correlation function to mea- 
sure H(z), Da (z), and fi(z). The obvious advantage of using mul- 
tipoles of the correlation function instead of the full 2D correlation 
function is the reduced number of data points used to obtain similar 
amount of information. 

The proper modeling of RSD is required in order to 
measure /3(z) or f B (z) from galaxy clustering data. Recent 
work on improving the modeli ng o f RSD include that of 
Ijennings. Baugh. & Pascolil ( 1201 ll) and IReid & White] fcOllh . In 
this paper, we focus on the detailed phenomenological modeling 
of the correlation function on smaller scales, in order to break the 
degeneracy between /3(z) and both H(z) and Da(z), to obtain sig- 
nificantly improved constraints on /3(z). We use the multipoles of 
the 2D correlation function for speed and efficiency. 

In Section[2] we introduce the galaxy sample used in our study. 
In Section[3] we describe the details of our new model. In Section 
[4] we test our model using mock data. In Section[5] we present our 
improved measurements from SDSS DR7 LRGs. We summarize 
and conclude in Sec. [6] 



2 DATA 

The SDSS has observed one-quarter of the entire sky and per- 
formed a redshift survey of galaxies, quasars an d stars in five 
passbands u, g, r, i, and z w ith a 2.5m telescope dFukugita et al.l 
ll996tlGunn et al.1 19981 l2006u . We use th e public catalog, the NYU 
Value-Added Galaxy Catalog (VAGC) dBlanton et aflfeoOSh , de- 
rived fr om the SDSS II final public data release, Data Release 7 
(DR7) ( Abazaii an et al. 2009). We select our LRG sample from 
the NYU VAGC with the flag primTarget bit mask set to 32. 
K-corrections have been applied to the galaxies with a fiducial 
model (ACDM with f2 m = 0.3 and h = 1), and the selected 
galaxies are required to have rest-frame g-band absol ute magni- 
tudes -23.2 < M g < -21.2 dBlanton & Roweisl2007l). The same 



selection criteria were used in previous paper s dZehavi et alj2005t 
lEisenstein et al.l 120051 ; lokumura et alj 120081 : IKazin et all l2010al) 
The sa mple we use is referred to as "DR7full" in IKazin et al.1 
d2010ij) . Our sample includes 87000 LRGs in the redshift range 
0.16-0.44. 

Spectra cannot be obtained for objects closer than 55 arcsec 
within a single spectroscopic tile due to the finite size of the fibers. 
To correct for these "collisions", the redshift of an object that failed 
to be measured would be assigned to be the same as the nearest 
successfully observed one. Both fiber collisi on corrections and K - 
corrections have been made in NYU- VAGC (Blan tonetaTll2005l) . 
The collision corrections applied here are different from what has 



been suggested in IZehavi et alj d2005l) . However, the effect should 
be small since we are using relatively large scale which are less 
affected by the collision corrections. 

We construct the radial selection function as a cubic spline fit 
to the observed number density histogram with the width Az = 
0.01. The NYU- VAGC provides the description of the geometry 
and completeness of the survey in terms of spherical polygons. We 
adopt it as the angular selection function of our sample. We drop 
the re gions with complet eness below 60% to avoid unobserved 
plates dZehavietai]|2005h . The Southern Galactic Cap region is 
also dropped. 



3 MODELING 2D CORRELATION FUNCTION 

In this section, we describe our model which encompasses the lin- 
ear scale to the nonlinear scale. 



3.1 Modeling 2D Correlation Function for Large Scales 

We comput e the linear matter power spectra , Pu n (k), by us- 
ing CAMB dLewis. Challinor. & Lasenby||2000l) . The linear power 



spectrum can be composed to two parts: 

Plin(k) =P nw (k) + P l B 



D Blo(k), 



(1) 



where P nw (k) is the no- wi ggle or pure CDM powe r spectrum cal- 
culated using Eq.(29) from lEisenstein & Hul d 19981) and PMao(^) 
is the wiggled part defined by the equation itself. The nonlinear 
damping effect o f the wiggled part in redshift sp ace can be well 
approximated by ( Eisenstei n. Seo. & Whit 



,(2) 



Isnitt spa 

ill! 



PBAo(k,p.k) = PsAo(k)-exp [ -^[1 + fi k (2f g + f g ) 



where fc» could be computed by dCrocce & Scoccimarrol 120061 
Matsubara 2008) 



1 

3^ 



1 -1/2 



Pu n (k)dk 



The dewiggled power spectrum is 

Pdiv(k,^k) = Pn-w(k) + PBAC>(k, /Ufc), 



(3) 



(4) 



fik is the cosine of the angle between k and the li ne of sight (LOS). 
Note t hat Eqs.([T}-(|4j are the same as Eq.(2) in IChuang & Wand 
d201 ll) . except for the addition of the direction-dependent terms in 
the exponent of the damping factor in Eq.(|2j, but are somewhat 
more intuitive. 

Next, we include the linear RSD as follows to obt ain the 
galax y power spectrum in redshift space at large scales d Raised 
L1983): 



Pg(k,fik) = b 2 (l + /3fil) 2 P dw (k,^i k ), 

= Pg,n-w(k,fJ.k) + Pg,BAo(k,^k), 

where b is the linear galaxy bias. Note that we have defined 

Pg,nw(k,Hk) = 6 2 (1 + /3nl) 2 P nw (k) 

Pg,BAo{k,ll k ) = b 2 (l + PfJ,k) 2 PBAo(k,lik) 



(5) 
(6) 

(7) 



= + p^t fPsAoik) ■ exp (^-^2 [1 + M*(2/ fl + fg])(8) 

Analogous to Eq.([6j(, the galaxy correlation function can be 
decomposed into no- wiggle and wiggled parts as follows: 
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€g,dw{<r,T) = ?g,n»(o",1") + £,"g,BAO 1") ■ 



(9) 



While Cg,d«j( (T ) 7r ) can De obtained by Fourier-transforming 
Pg (k, fik), doing so involves two-dimensional integrals, and thus is 
time-consuming and inefficient. Instead, we can Fourier transform 
each term in Eq.JfjJ separately, using Legendre polynomial expan- 
sions and integral convolutions that only involve one-dimensional 
integrals. 

The no-wiggle galaxy correlation function in redshift space 
can be compute d by Fourier transforming Eq.lO, which gives 
dHamiltoiJI 19921) 

= b 2 (ed W is)P (fi)+Q w (s)P2(fi)+a w is)P i (fi)),m 



where s = \/a 2 + tv 2 , li is the cosine of the angle between 
s = (a, 7t) and the LOS, and Pi are Legendre polynomials. The 
multipoles of are defined as 



Co (r) 
2 (r) 
U (r) 



35 



7— nit; 

|C (r) 



where /? is the linear RSD parameter and 



(r) = 

— niw 

£ M = 



/-rati; / l\ ll 7 / 

£ (r )r ar , 

j-TLW f l\ I A j I 



(ii) 

(12) 
(13) 

(14) 
(15) 



where £ n ™(r) is obtained by fourier transforming Pnw(k). 

The wiggled part of the galaxy correlation function in redshift 
space, Cg,sAo( cr ! 7r )' is obtained by Fourier transforming Eq.©. 
Note that the ^t-dependent damping factor in fc-sapce in Eq.([8jl be- 
comes a Gaussian convolution in configuration space: 



Cg,BAo{ a : ^) 



£*(<7, TV — x) f*(x)dx, 



(16) 



where (a, ir) is the Fourier transform of b 2 (1 + /3liI) 2 Pbao ' 
cxp(-^-), and 

/*(z) = _!_exp(— V (17) 

a-ty/TY \ at 1 

where 



2 _ 4/ g + 2/ 2 



(18) 



£* (cr, 7r) can be obtained using Eq. 4 1 01>-d 15b. but replace £ m0 (r) 
(the Fourier transform of P nU ;(fc)) with the Fourier transform of 

P»2 (fe)-exp(^). 



galaxy bias at small scales by multiplying £ nw (r) with the follow- 
ing factor 



(19) 



where 6 a is a constant. F(r) is a function which is close to 1 for 
small r and close to when r is large; we choose 



F(r) = 



(20) 



where we choose 6s = 30/i -1 Mpc and &c = 4; these are 
motivated by the fact that the galaxy correlation function is a 
powerlaw at small scales (i.e. s < 15/i _1 Mpc) and the scale- 
dependent effects (including nonlinear effects and scale-dependent 
galaxy bias) are negligible at larger scales, s > 40/i~ 1 Mpc. The 
overall scale-dependent effects are included when computing the 
no-wiggle galaxy correlation function by replacing £ nw (r) with 
£ nw (r) x b n i(r) in applying Eq. lfT0t-lfT5t. The resultant corre- 
lation function is denoted as £,g'^ l w (a, ir). 

We now obtain the 2D correlation function that incorporate 
nonlinear effects, galaxy bias, and linear RSD: 



(21) 



where £g i BAa{ a i 7r ) K given by Eq.lfTBl. 

Next, we convolve the 2D correlation function with the distri- 
bution function of r andom pairwi se velocities, f(v), to obtain the 
final model f (<r, ir) dPeebleslll980l) 



£0, ir)= f cr, 7T 



ff(»a(» 



f(v)dv, 



(22) 



where the random motions are rep resented by an exponential form 
jRatcliffeet all 19981 ; lLandvf2 002) 



/(«) 



1 ( y/2\v\ 

■. exp ' 



(23) 



where a v is the pairwise peculiar velocity dispersion, s ■ = a + 

{ir- 



« ^2 „„ . ,./ _ 1 



and jix = -y(7r ■ 



,, — , ,„ u, w \ J ■ We find that the 2D 
correlation functions measured from LasDamas mocks can be well 
fitted by 



a v (s',Li' 2 ) = <T„, (1 + c MA t ,2 )(l + c CT e <nsr), 



(24) 



where ovo is the dispersion corresponding to the truly random mo- 
tion and c M and c a terms describe the dependence on direction and 
se paration respectively. The q -dependence is similar to that found 
bv lCabre & Gazta naea ( 2009). They found that the 2D correlation 
functions from the MICE N-body simulations are fitted well with a 
pairwise velocity distribution which is large when a < 5/i _1 Mpc. 
We have added the direction-dependent term, c M /i' 2 , to model the 
high amplitude of £4 at small scales (see Fig.|4j- 



3.2 Modeling 2D Correlation Function for Small Scales 

For small scales, we need to model three effects: the nonlinear mat- 
ter correlation function, the scale-dependent galaxy bias, and the 
RSD from the random galaxy pairwise velocities. It is well known 
that the small sc ale galaxy co rrelation function is well described 
by a powerlaw aPeeblesll 19801) . Since the galaxy correlation func- 
tion is given by £ nw (r) on small scales, we model the combination 
of nonlinear matter correlation function and the scale-dependent 



4 VALIDATION OF THE MODEL USING MOCK 
CATALOGS 

In this section, we present the methodology and results of testing 
our model described in the previous section. 
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4.1 Mock Catalogs Used 

We use the 160 mock catalogs from the LasDamas simulation^ 
(McBride et al., in preparation) to test our model. LasDamas pro- 
vides mock catalogs matching SDSS main galaxy and LRG sam- 
ples. We use the LRG mock catalogs from the LasDamas gamma 
release with the same cuts as the SDSS LRG DR7full sample, 
-23.2 < M 3 < -21.2 and 0.16 < z < 0.44. We have diluted the 
mock catalogs to match the radial selection function of the obser- 
vational data by randomly selecting the mock galaxies according to 
the number density of the data sample. We calculate the multipoles 
of the correlation func tions of the mock catalo gs and construct the 
covariance matrix (see lChuang & Wangl j2012h for details). 



4.2 Measuring the Two-Dimensional Two-Point Correlation 
Function 



!. ,'■ v. . *. 

>i 




-150 -100 -50 50 100 150 

Transverse Separation (Mpc/h) 



We convert the measured redshifts of galaxies to comoving dis- 
tances by assuming a fiducial model, ACDM with Q m — 0.25. 
We use the two-poin t correlation function estimator given by 
lLandv & SzaiavUl993l) : 



4>,tt) = 



DD(a, tt) - 2DR(a, it) + RR{o, ir) 
RR(a, tt) 



(25) 



where n is the separation along the line of sight (LOS), a is the sep- 
aration in the plane of the sky, DD, DR, and RR represent the nor- 
malized data-data, data-random, and random-random pair counts 
respectively in a given distance range. The LOS is defined as the di- 
rection from the observer to the center of a pair. The bin size we use 
here is 1 fa _1 Mpcx 1 h~ 1 Mpc. The Landy and Szalay estimator has 
minimal variance for a Poisson process. Random data are generated 
with the same radial and angular selection functions as the real data. 
One can reduce the shot noise due to random data by increasing the 
number of random data. The number of random data we use is 10 
times that of the real data. While calculating the pair counts, we as- 
sign to each data point a radial weight of 1/[1 + n(z) ■ P w ], where 
n(z) is the radial sele ction function and P w = 4 ■ 10 4 h~ 3 Mpc' i 
jEisenstein et alj|2005h . 

Fig [T] shows the averaged 2D correlation function measured 
from the mock catalogs. We use the averaged radial selection func- 
tion to construct the random catalog since it is closer to the true 
mean density. Clearly, our model provides an excellent fit to data 
over a wide range of scales, from the largest scales where data are 
not too noisy, to the smallest scales plotted (except very near the 
LOS). 



4.3 Multiples of the Correlation Function 

As in IChuang & Wangl J2012h . the effective multipoles of the cor- 
relation function are defined by 



(2/ + l)C(<7,7T)P i ( M )Vl-M 2 



^ /g\ (26) 

Number of bins used in the numerator 

where As = 5 /i _1 Mpc in this work, and 

a= (n+ i)ft _1 Mpc,n = 0, 1,2, ... (27) 
7r = (m+i)/i" 1 Mpc I m = 0,1,2,... (28) 



Figure 1. The average two-dimensional two-point correlation function (2D 
2PCF) measured from 160 LasDamas SDSS LRGfull mock catalogs (solid 
black contours), compared to a theoretical model with the input parameters 
of the LasDamas simulations and 64, er„ n, Cfi, Co-} are set to {0.316, 
—0.0385, 50km s —1 , 10, 4} (dashed red contours). The thick dashed blue 
circle denotes the baryon acoustic oscillation scale. The contour levels are 
£ = 2,0.5,0.1,0.025,0.01,0.005,0. The § = contours are denoted 
with dotted lines for clarity. 



/'• : - 



y<7 2 + n 2 • 



(29) 



Note that both the measurements and the theoretical predic- 
tions for the effective multipoles are computed using Eq.J26t. We 
do not use the conventional definitions of multipoles to extract pa- 
rameter constraints as they use continuous integrals. Bias could be 
introduced if the definitions of multipoles are different between 
measurements from data and the theoretical model. 

Fig-IUm and|4]show the effective monopole (£0), quadrupole 
(£2), and hexadecapole (£4) measured from the LasDamas mock 
catalogs comparing to our full model and a simpler model (linear 
model + ID dewiggle damping + constant velocity dispersion). In 
Fig. [2] one can see how our model completely corrects the scale- 
dependent effects in the measured monopole. Fig[3] shows that our 
model provides a reasonable fit to the measured quadrupole. In Fig. 
[4] we find that angle-dependent term, c CT , significantly improves the 
fitting of hexadecapole at small scales (s < 50/i _1 Mpc). However, 
at larger scales(s > 60h~ Mpc), the LasDamas mocks show some 
oscillatory features while the theoretical models are flat. It is likely 
due to the dewiggle damping not being adequate enough to model 
£4 and one might need higher order term (i.e. /1 4 ). Therefore, we 
do not include £4 to measure parameters in this study. 



4.4 Covariance Matrix 



We construct the covariance matrix as 



A 



(30) 



where N is the number of the mock catalogs, X m is the mean of 
the m th element of the vector from the mock catalog multipoles, 
and X^ is the value in the m th elements of the vector from the k th 
mock catalog multipoles. The data vector X is defined by 



http://lss.phy.vanderbilt.edu/lasdamas/ 



(1) tf(2) 



i(JV).£(l) i(2) 
l SO ' S2 ' S2 i 



•I S2 1 



(31) 
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Figure 2. The averaged monopole of the correlation functions of the mock 
catalogs (black squares) comparing to the fitting model of this study (red 
dots) and a simpler model (linear model + ID dewiggle damping + constant 
velocity dispersion, blue crosses). The error bars are taken as l/\/160 of 
the square roots of the diagonal elements of the covariance matrix. 



Figure 4. The averaged hexadecapole of the correlation functions of the 
mock catalogs (black squares) comparing to the fitting model of this study 
(red dots) and a simpler model (linear model + ID dewiggle damping + con- 
stant velocity dispersion,blue crosses). The error bars are taken as 1/ \J 160 
of the square roots of the diagonal elements of the covariance matrix. 



LasDamas Mocks 
Fitting model in this study 

Lineal" model + ID dewiggle + constant velocity dispersion 



where Nx is the length of the vector used, X t h is the vector from 
the theoretical model, and X t s is the vector from the observational 
data (we use the mock catalogs as the observational data to test the 

model in this section); 

As explained in Chu ang & Wan3 ( l201 ll) . instead of recalculat- 
ing the observed correlation function for different theoretical mod- 
els, we rescale the theoretical correlation function to avoid render- 
ing x 2 values arbitrary. The rescaled theoretical correlation func- 
tion is computed by 



It 



V. 



Tt ' fT inT = f H^''-ffW" 



H(z) 

where S^th is given by eq. i22l . Hence \ 2 can t> e rewritten as 



(33) 



20 40 60 80 100 120 

s (Mpc/h) 



Figure 3. The averaged quadrupole of the correlation functions of the mock 
catalogs (black squares) comparing to the fitting model of this study (red 
dots) and a simpler model (linear model + ID dewiggle damping + constant 
velocity dispersion, blue crosses). The error bars are taken as 1/V160 of 
the square roots of the diagonal elements of the covariance matrix. 



where N is the number of data points in each measured multipole; 
TV = 19 while using the scale range, 25 < s < 120/i _1 Mpc. The 
length of the data vector X depends on how many multipoles are 
used. 



4.5 Likelihood 

The likelihood is taken to be proportional to exp(— x 2 /2) 
jPress et alii 19921) , with x' 2 given by 



N X 

E 



yfid 1 
oba.j { 



(34) 



^ ' [Xth,i ~ X bs,i] Ctj 1 \Xth,j — X b a ,j\ 



(32) 



where T~ 1 X t h is a vector given by eq. d 3 3 1 > with £ t h replaced 
by its effective multipoles (defined by eq. l!26ll), and X^f is the 
corresponding vector from observational data measured assum- 
ing the fiducial model in converting redshifts to distances. See 
IChuang & Wan3 j201ll) for a more detailed description of our 
rescaling method. 



4.6 Markov Chain Monte-Carlo Likelihood Analysis 

We use CosmoMC in a Mark ov Chain Monte-Carlo likelihood 
analysis dLewis & Bridie1l2002l) . The parameter space that we ex- 
plore spans the parameter set of {i7(0.35), Da(0.35), Q m h 2 , 
13, Q b h 2 , ris, b A , (7 v ,o, c M , c CT }. Only {//(0.35), Da(0.35), 
Q m h 2 , /?} are well constrained by the mock data. We marginal- 
ize over {6a, <?v,o, c m , c ct } with the flat priors, {(—0.2,0.2), 
(0, 500)s _1 km, (0,20), (0,10)}. We marginalize the values of 
Q.bh 2 and n s with a wide range (±7ctwmap iKomatsu et al.l 
(2010)) centered at the input values of simulations (for mocks) or 
WMAP7 measurements. We also marginalize over the amplitude of 
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the galaxy correlation function. For simplicity, we fix f g (z = 0.35) 
in the asymmetric damping factor in Eq.l[2} to the fiducial value, 
0.645. Our results are not sensitive to this choice. 

4.7 Validation of the Model Using One Mock 

To simulate the application of our model to the SDSS DR7 LRGs, 
we apply it to one LasDamas SDSS LRG mock catalog, to test 
whether the true values of the parameters (the input parame- 
ters) are recovered by our methodology. Table Q] shows the mea- 
surements of {#(0.35), £>a(0.35), n m h 2 , /3, #(0.35) r 3 {z d )/c, 
Da(0.35) /r s (z d )} from one LasDamas SDSS LRG mock catalog 
(indexed Ola in the LasDamas simulations). Clearly, when smaller 
scales are included, the constraints are improved, while the mean 
values generally become closer to the input values (although there 
are fluctuations due to sample variance). Finally, we find that x 2 
per degree of freedom is reasonable, and reflects the level of noise 
and systematic bias inherent in the mock data. 



5 IMPROVED MEASUREMENTS FROM SDSS DR7 LRG 

Table [2] lists the mean and rms variance of the param- 
eters, {i/(0.35), £>a(0.35), Q m h 2 , /3, #(0.35) r s (z d )/c, 
Da(0.35) /r s (z d )}, derived in anMCMC likelihood analysis from 
the measured £o + £2 of the correlation function of the SDSS LRG 
sample with the scale range, 25 < s < 120 /i _1 Mpc. Table|3]gives 
the corresponding normalized covariance matrix. 

We do not use smaller scales (i.e. 15 < s < 120/i _1 Mpc) 
because it gives a much too high x 2 /d.o.f. (larger than 2). We 
marginalize over the other parameters, {Qth 2 , n s , cr v }, with flat 
priors over the ranges of {(0.01859,0.02657), (0.865,1.059), 
(0, 500)s~ 1 km}, where the flat priors of Qbh 2 and n a are cen- 
tered on the measurements fro m WMAP7 and has width of 
in-law map (with aw map from iKomatsu et aL These 
priors are wide enough to ensure that CMB constraints are not 
double counted when our results a re combined with CMB data 
dChuang. Wang. & Hemantha 2012). We also marginalize over the 
amplitude of the g alaxy correlation functio n. 

Our model in lChuang & Wand J2012I) was only applicable for 
40 < s < 120fo -1 Mpc. Comparing our results from this study 
(25 < s < 120fc~ M pc) and those from our previous study in 
IChuang & Wmgj J2012h (40 < s < 12Q/i _1 Mpc) in Tabled one 
can see the dramatic improvement in the constraints obtained from 
this study from adding the information on smaller scales. 

It is important to note that the correlation between 
£(0.35) and {#(0.35)r s (z d )/c, D A (0.35)/r s (z d )} are 
significantly reduced: r(j9(0.35), H(Q.35)r s (z d )) = 
0.748 and r(/3(0.35), D A (0-35)/r s (z d )) = 0.575 fo r 
40 < s < 12Q/i _1 Mpc l lChuang & Wmgj |2012|) , 
versus r(/3(0.35), #(0.35)r s (z d )) = 0.494 and 
r(/3(0.35),D A (0.35)/r s (2 d )) = 0.159 for 25 < s < 
120/i _1 Mpc (this study). 





25 < 8 < 120 


40 < s < 120 


ff(0.35) 


76.3 ± 5.6 


79.6 ±8.8 


D A (0.35) 


1126 ± 49 


1060 ± 92 


1 Ifji il 


0.097 ± 0.012 


0.103 ± 0.015 





0.348 ± 0.058 


0.44 ± 0.15 


H(0.35)r 3 (z d )/c 


0.0423 ± 0.0027 


0.0435 ± 0.0045 


D A (0.35)/r s ( 2d ) 


6.77 ± 0.17 


6.44 ± 0.51 


X 2 /d.o.f. 


1.65 


1.23 



Table 2. The mean and standard deviation of {i7(0.35), D A (0.35), 
tt m h 2 , /3, H(0.35)r 3 {z d )/c, D A (0.35)/r s (z d )} from SDSS DR7 
LRGs using £0 + £2 and the scale range, 25 < s < 12 0fe~ 1 Mpc, com- 
parin g with the measurements from our previous paper jChuang & Wand 
120121) using only large scales, 40 < s < 120/i~ 1 Mpc. The x 2 /d.o.f. is 
1 .65 which is a bit high but still could be due to the sample variance. The 
unit of H is km s — 1 Mpc -1 . The unit of Da and r B (z d ) is Mpc. 

Applying this model to the SDSS LRGs over the scale range of 
25 < s < 120/i _1 Mpc, we are able to break the degeneracy be- 
tween f3(z) with both H(z) and Da(z), and obtain significantly 
improved measurements: H(z)r s (z d )/c = 0.0423 ± 0.0027, 
D A {z)/r s (z d ) = 6.77 ± 0.17, and /3(z) = 0.348 ± 0.058 at 
z = 0.35. 

Our model incorporates the overall nonlinear effects via 
the use of the "dew iggled" galaxy power spectrum, as in 
IChuang& Wand" d20 111) , but we now include the enhanced damp- 
ing along the line of sight (see Eqs.lfT)-©), even though it only 
has a small effect on the results. On small scales, the nonlinear ef- 
fect and scale-dependent galaxy bias are degenerate, and we model 
these as an overall scale-dependent correction. Most significantly, 
we allow the RSD to be scale and direction dependent in our model. 
Our model provides excellent fit to mock data (see Fig[T]l. 

In order to test our model in a timely fashion, we have cho- 
sen to use only the monopole and quadrupole of the 2D correlation 
function, instead of the full 2D correlation function, in our anal- 
ysis of the mock and real data. While our model provides an ex- 
cellent fit to mock data down to a scale of 15 /i -1 Mpc, it leads to 
a steep increase in the x 2 per degree of freedom at scales smaller 
than 25 /t^Mpc when applied to the SDSS DR7 LRGs. However, 
the constraints on H(z)r s (z d ), DA(z)/r s (z d ), and j3{z) remain 
consistent at the 68% C.L., while the errors on these increase as 
we include scales below 25 /i _1 Mpc. This may indicate that the 
full 2D correlation function (instead of just its monopole and the 
quadrupole) needs to be used to model the data on smaller scales. 
We will do so in our future work. We will apply our model to the 
public BOSS data, which will allow a greater range of scales (on 
both small and large scales) to be used in the analysis for measuring 
H(z)r s (z d ), D A {z)/r s (z d ), and f3{z). 

We expect our model to be useful in tightening dark energy 
and gravity constraints from the full analysis of current and future 
galaxy clustering data. 



6 CONCLUSION AND DISCUSSION 

We have presented and validated a simple and efficient phenomeno- 
logical model for the two-dimensional two-point galaxy correla- 
tion function that works well over a wide range of scales, from 
large scales down to small scales not used in our previous work 
(where we restricted ourselves to scales larger than 40/i~ 1 Mpc). 
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15 < s < 120 


25 < s < 120 


30 < s < 120 


40 < s < 120 


input value 


ff(0.35) 


86 ± 11 


78 ± 13 


79 ± 12 


83 ± 13 


81.79 


Da (0.35) 


1050 ± 92 


1003 ± 112 


1008 ± 110 


1035 ± 116 


1032.8 


O U 2 


0.125 ± 0.018 


0.115 ± 0.017 


0.114 ± 0.018 


0.116 ± 0.019 


0.1225 





0.41 ± 0.12 


0.44 ± 0.15 


0.50 ±0.16 


0.57 ± 0.17 




H (0.35) r s (z d )/c 


0.0457 ± 0.0052 


0.0423 ± 0.0064 


0.0429 ± 0.0062 


0.0451 ± 0.0068 


0.0434 


D A (0.35)/r s (z d ) 


6.59 ±0.54 


6.17 ± 0.67 


6.18 ±0.64 


6.36 ± 0.68 


6.48 


X 2 /d.o.f. 


1.24 


1.09 


1.06 


1.03 





Table 1. The mean and standard deviation of the distribution of the measured values of {H(0.35), Da (0.35), tt m h 2 , 0, H (0. 35) r a (z d ) / c, 
D A (0.35) /r a (z d )} from one LasDamas SDSS LRG mock catalog (indexed Ola) using four scale ranges, 15 < s < I20h~ Mpc, 25 < s < 120/i -1 Mpc, 
30 < s < 120/i -1 Mpc and 40 < s < 120fe~ 1 Mpc. The constraint of is significantly improved when including smaller scales. The mean value of is 
also closer to expected value, ~ 0.32 (see lChuang & W ang 1 2012)). The unit of the separation, s, is h — x Mpc. The unit of H is km s — 1 Mpc - 1 . The unit of 
D A and r a (z d ) is Mpc. 





H(0.35) 


Da (0.35) 


m h 





H (0.35) r a (z d )/c 


D A (0.35)/r a (z d ) 


H(0.35) 


1.0000 


-0.2788 


0.4525 


0.4160 


0.8625 


0.2597 


D A (0.35) 


-0.2788 


1.0000 


-0.7617 


0.1112 


0.1515 


0.5092 




0.4525 


-0.7617 


1.0000 


-0.0856 


-0.0290 


0.0869 





0.4160 


0.1112 


-0.0856 


1.0000 


0.4941 


0.1594 


H(0.35)r s (z d )/c 


0.8625 


0.1515 


-0.0290 


0.4941 


1.0000 


0.2483 


D A (0.35)/r 3 (z d ) 


0.2597 


0.5092 


0.0869 


0.1594 


0.2483 


1.0000 



Table 3. Normalized covariance matrix of the measured and derived parameters, {H(0.35), Da (0.35), O m h 2 , 0, H(0.35) r s (z d )/c, D A (0.35) /r a (z d )} 
from SDSS DR7 LRGs using | + |a and the scale range, 25 < s < 120/i -1 Mpc. 
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